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The coupling between mechanical and thermal properties due to thermal expansion complicates 
the problem of measuring frequency-dependent thermoviscoelastic properties, in particular for highly 
viscous liquids. A simplification arises if there is spherical symmetry where - as detailed in the 
present paper - the thermoviscoelastic problem may be solved analytically in the inertia- free limit, 
i.e., the limit where the sample is much smaller than the wavelength of sound waves at the frequencies 
of interest. As for the one-dimensional thermoviscoelastic problem [Christensen et al., Phys. Rev. E 
75, 041502 (2007)], the solution is conveniently formulated in terms of the so-called transfer matrix, 
which directly links to the boundary conditions that can be experimentally controlled. Once the 
transfer matrix has been calculated, it is fairly easy to deduce the equations describing various 
experimentally relevant special cases (boundary conditions that are adiabatic, isothermal, isochoric, 
etc.). In most situations the relevant frequency-dependent specific heat is the longitudinal specific 
heat, a quantity that is in between the isochoric and isobaric frequency-dependent specific heats. 

PACS numbers: 64.70.P- 



I. INTRODUCTION 

Linear thermoviscoelasticity is the well-established discipline dealing with the irreversible thermodynamics of slightly 
perturbed systems where mechanical and thermodynamic degrees of freedom couple to each other d, 01 • For glass- 
forming liquids cooled towards the calorimetric glass transition relaxation times become very long compared to phonon 
times, approaching and eventually exceeding seconds. For such "uitraviscous" liquids not only the mechanical moduli 
become complex and frequency dependent [1, @, Q) but so do standard thermodynamic linear-response properties like 
i the specific heat or the thermal expansion coefficient [1, i, [M HH, E E III, [ll [ll, [11, [ll H [lE [12, [B, |M 111, 
' Of the twelve basic c omp lex, frequency-dependent thermodynamic linear-response coefficients only three are inde- 
pendent (see, e.g., Ref. [23 and its references). If one assumes stochastic dynamics, which is believed to be realistic for 
viscous liquids on time scales much longer than phonon times, in fact only two thermoviscoelastic response functions 
are truly independent [l], H, [13, [II]- No reliable measurements of a full set of (three) thermoviscoelastic response 
function appear to exist for any highly viscous liquid. Part of the reason for this may be the traditional focus in 
physics on phenomena on the molecular scale, but part of the problem most likely also comes from the considerable 
challenges associated with reliably measuring the thermoviscoelastic response functions. The problem is that the 
' coupling between mechanics and thermodynamics caused by thermal expansion is nontrivial when the mechanical 
shear modulus is a significant fraction of the bulk modulus [1^. This is true for solids in general, as well as for 
highly viscous liquids at frequencies of order the inverse relaxation time. For solids, however, because the isobaric 
and isochoric specific heats are almost identical, the problem is not serious. For uitraviscous liquids, on the other 
hand, the problem cannot be ignored. For such systems it was recently shown that the thermomechanical coupling 
implies that conventional methods fail to measure the isobaric, frequency-dependent specific heat 23]. This is because 
truly isobaric conditions are difficult to establish experimentally due to the thermal expansion upon heating; in most 
I experiments attempting to measure the frequency-dependent specific heat the stress tensor is not proportional to the 
unit tensor (i.e., there is not hydrostatic equilibrium), and shear stresses relax on the very time scale that one wishes 
to monitor. 

The purpose of this paper is to establish the theoretical framework for measuring a complete set of thermoviscoelastic 
response functions utilizing spherical symmetry. The thermomechanical equations describing the coupling of mechanics 
to thermodynamics are well known since many years [l[,[3. We recently (with Olsen) presented the full analytic solution 
of the one-dimensional inertia- free case, i.e., where all motion is restricted to one direction and the wavelength of sound 
waves at the relevant frequencies is much larger than the sample size j23j . The reader is referred to this paper as an 
introduction to the reasoning and techniques used in the present paper. Below, the spherically symmetric inertia-free 
case is treated, where the effective one-dimensional nature again allows the problem to be solved analytically. The 
solution is cumbersome, but once it has been arrived at a number of experimentally relevant special cases are fairly 
easy to calculate. 

The solution is formulated in terms of the so-called transfer matrix (23l . [soj that links infinitesimal variations at 
the boundaries for the following four quantities: entropy (or equivalently heat) input, volume displacement, pressure 
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FIG. 1: Energy bond graphs [3ll . [3^ |33| . [sj], a useful tool for modeling linear thermoviscoelasticity. (a) The standard 
thermodynamic energy bonds, one thermal and one mechanical. In energy bond terminology temperature and negative pressure, 
T and —p, are so-called effort variables analogous to voltage, and the entropy flux S = dS/dt and volume flux V = dV/dt are 
so-called flow variables analogous to electrical current. In quasi-equilibrium the product of an effort variable and its flow variable 
is the energy flux into the system from its surroundings; more generally it gives the flux of free energy - in the sense of available 
work - into the system, (b) Symbolic flgure of the energy bonds of linear irreversible thermodynamics for the spherically 
symmetric situation. For the thermal energy bonds the efforts are the temperature variations STi and 5T2 at radius ri and 
r2, respectively, and the flows are the entropy fluxes SSi and -682 into the system (fluxes are in the positive radial direction). 
For the mechanical energy bonds the flows are the volume fluxes —SVi and SV2 into the system, respectively, whereas pressure 
is now the "radial pressure," Spr,i and 5pr,2 respectively, deflned in Eq. (|44)l . The four energy bonds deflne eight variables. 
As shown in this paper, the fundamental physical equations provide four constraints among these eight variables. Thus there 
is a linear relation between the four "outer" and the four "inner" energy bond variables. This relationship is expressed in 
terms of the 4x4 transfer matrix calculated below that provides all information needed to interpret any experimental situation 
characterized by specific boundary conditions. 

change, and temperature change. The transfer matrix - not to be confused with the transfer matrix of statistical 
mechanical models - is useful because it directly describes how the system interacts with its surroundings. These 
interactions are conveniently pictured and described via the energy bond technique [sil [s^ . [ssl . [34l] . An energy bond 
is characterized by a displacement variable and an effort variable, generalizing the concepts of charge and voltage. 
The product of an effort and a differential displacement variable is a generalized work that gives the (free) energy 
transferred into the system from its surroundings. 

In standard thermodynamics there are two energy bonds, one thermal and (if the shear modulus is negligible) one 
mechanical, see Fig. 1(a). The thermal energy bond is characterized by entropy S as the displacement variable and 
temperature T as the effort variable. For the mechanical energy bond the displacement is the volume V ^ and the 
effort is the negative pressure, —p. The time derivatives of displacement variables define flow variables, signaled by 
arrows in Fig. 1. The two energy bonds of Fig. 1(a) represent the fundamental thermodynamic identity 



In linear irreversible thermodynamics small perturbations in temperature and pressure are linear functionals of the 
volume and entropy flows, or vice versa. The temperature perturbation dT = T — Tq is around a reference temperature 
To, and the pressure perturbation dp = p — pq around a reference pressure pq. The supplied heat dQ is related to the 
externally supplied entropy, dScxt = dQ/T. When the volume Vq is so small that temperature and pressure can be 
assumed homogenous throughout the volume, the relaxation is described in the time domain by memory kernels of 
the form (where kt is isothermal compressibility, ap isobaric expansion coefhcient, and Cp isobaric specific heat per 
unit volume) 



dE = TdS - pdV . 



(1) 
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dV{t)lVQ = - / nT{t-t')dp{t')+ j ap{t-t')dT{t'), (2) 

'OO 'J —OO 

t rt ^ 



dS,^t{t)/Vo = - ap{t - t')dp{t') + Cp{t-t')dT{t'). (3) 

J-oo J-oo -'0 

Instead of Eq. ([J) we now have 

dE ~ TdS'oxt — pdV — dTdSoxt — dpdV + TodScxt ~ PodV. (4) 

or 

d{E ~ ToS-cxt + PoV) = dTdScxt - dpdV. (5) 

The function E — ToScxt + PoV is the imparted free energy, which is not a state function since dissipation in the 
system degrades this energy: During a cychc process one has 



diE - To5ext + PqV) = -To i dS,xt = TqS , (6) 



where S is the entropy production in the system. Note that the entropy production is quadratic in the perturbations, 
implying that entropy is conserved to first order a fact that is utiUzed below. 

It follows from the above that the generic conjugated variables for a small system that relaxes mechanically and 
thermally with no shear forces are (dT, diSext) and {—dp, dV). From an experimental point of view, however, it is not 
convenient only to consider an infinitesimal volume, Vq. If one wishes to study relaxation on a time scale Toxp, the 
volume Vq may only be considered small if the heat diffusion time td across the volume is much smaller than Texp- 
This is an unduly restriction that can only be coped with by analyzing the influence of heat diffusion on the response - 
an important purpose of this paper. A further complication arises when the shear modulus becomes comparable to the 
bulk modulus, which is the case in the relaxation region of viscous liquids. In this case heat diffusion and mechanical 
stresses couple nontrivially. To keep the below discussion as simple as possible we look at the situation with highest 
symmetry, that of a sphere of inner radius ri and outer radius r2, posing the question: What is the relation between 
the thermal and mechanical variables at the boundaries of the system? In order to take the shear forces properly into 
account it is shown below that pressure must be replaced by the "radial pressure," Spr,i at ri and (5pr,2 at r2- 6Vi 
is the volume swept by the surface at ri as a consequence of a radial, small displacement u{ri). Correspondingly, 
SV2 is the volume swept at r2 (both in the positive radial direction). The net volume change is dV = —SVi + 6V2- 
The temperature perturbations at the two surfaces are denoted by STi and ST2, respectively, and the entropy fluxes 
in positive radial direction by 6 Si and SS2, respectively. The net entropy influx is dS'ext = — 682- Now Eq. ^ 
becomes 



d{E - ToS,xt+PoV) = STiSSi - 5T25S2 + SprsSVi - 5pr,2SV2. 



(7) 



As illustrated in Fig. 1(b) this gives rise to four energy bonds, two referring to the outer radius r2 and two to the 
inner radius ri. 

The transfer matrix T(r2,ri) is by definition the 4x4 matrix that links the (generally complex and frequency- 
dependent, see below) infinitesimal variations of the four energy bond variables at radius ri with the four energy bond 
variables at radius r2. Switching notation such that 6pr 1 is denoted by Sp{ri), etc, the transfer matrix is thus defined 
by 



/ 5p{r2)\ / 5p[ri) 

5T{r2) A ST{n) 

\SSir2)/ \SS{n), 



(8) 



The transfer matrix is related to the response matrix r(r2,ri) that by definition links the four effort variables to the 
four displacement variables: 
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ST{n) 
5p{r2) 



= r(r2,ri] 



/ SV{n) 
SSin) 
-SV{r2) 



(9) 



From the fluctuation-dissipation theorem the response matrix is known to be symmetric, a fact that is exphcitly 
confirmed below. 

In the next section fundamentals are summarized. In Sec. Ill the full dynamic equations are formulated and 
brought into dimensionless form by scaling with complex units. In Sec. IV the equations are solved, and the transfer 
and response matrices are calculated. In Sec. V several experimentally relevant special cases are considered. Finally, 
Sec. VI gives a brief discussion. 



II. DEFINITIONS AND CONSTITUTIVE RELATIONS 

The energy bond variables of standard thermodynamics give rise to a number of dc (i.e., static) linear-response 
coefficients as follows. If the variables of interest are those of the two thermodynamic energy bonds of Fig. 1(a), 
{T ,p ,S , V), there are altogether 24 thermodynamic coefficients of the form {da/db)c with a, b, and c chosen among 
T,p, S, and V [25l[35j. These coefficients form 12 pairs that are trivially related by inversion [{da/db)c = l/{db/da)c, 
etc]. As is well known the 12 coefficients are not independent, but related by Maxwell relations. This leaves the 
following eight basic linear-response coefficients (where the specific heats here and throughout the paper are per unit 
volume): 

Isochoric specific heat: cy = (10) 

Isobaric specific heat: Cp = -^f^^) (11) 

1 /'dV\ 

Isothermal compressibility: kt ^ ( (12) 

V \ dp / T 

1 /-dVx 

Adiabatic compressibility: ks = ~Y\~Q~)g ^^'^^ 

T , • • rr^ ■ 1 fdy\ 1 /dS\ , 

Isobaric expansion coefficient: a„ = — = t;— (14) 

V \ dT ) p V \dp/T 

, , • • ■ 1 fdV\ 1 fdS\ 

Adiabatic contraction coefficient: as = ^ — = — -77- 15) 

V\dTJs V\dpJv ^ ' 



Isochoric pressure coefficient: By = f = 

\oT/v \dV/T 

Adiabatic pressure coefficient: 63 = { tt- ) = ( — 1 

\dT/s \dV J p 



(16) 
(17) 



Some well-known relations between the thermodynamic coefficients are summarized in the Appendix for reference, 
where relations are simplified somewhat by letting the heat capacities be represented by the related variables Qv = 
1 C^'S'It, and r = ^ 

In systems with relaxing degrees of freedom the thermodynamic coefficients generally become complex and frequency 
dependent. Suppose, for instance, that the system is subjected to an infinitesimal periodic pressure variation with 
angular frequency w described as p{t) = pa + Re[(5pexp(si)], where s = ±.iuj depending on convention [g^l is the so- 
called Laplace frequency. The volume then varies periodically as V{t) = Vq + Re[(5Fexp(si)], where both 5p and 5V 
are generally complex and frequency dependent. If this takes place at constant temperature, the complex frequency- 
dependent isothermal compressibility is defined via kt = —5V/ (JpVo), where Vq is the average volume. According to a 
basic theorem of linear irreversible thermodynamics, the Maxwell relations among the dc linear-response quantities in 
Eqs. (6)-(9) translate into Onsager relations for the frequency-dependent coefficients, which reflect time reversibility. 
This is a special case of the so-called correspondence principle [U, [2, @| ^ Any (dc) thermodynamic relation or equation 
involving linear thermodynamic and/or mechanical quantities applies unchanged when constitutive properties are 
replaced by the corresponding complex, frequency-dependent quantities. These frequency domain functions are related 
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to the corresponding memory kernels like, e.g., 

j'OO 

kt{s) =s HT{t)e^'*dt, (18) 

^0 

where we adhere to the ordinary (sloppy) notation of physics that it is to be read out of the context whether kt refers 
to time or frequency domain function. In the following we stay in the frequency domain, though. 

Thermoviscoelasticity describes the coupling between thermal and mechanical deviations from equilibrium. This 
paper deals with the linear case that is well understood as regards fundamentals. Linearity means that the system is 
assumed to be infinitesimally close to equilibrium. Deviations from equilibrium are quantified in terms of the infinites- 
imal displacement field u = u(r, t), temperature variation field ^T(r, t) = T{r, t) — Tq, etc. In this approximation the 
radial heat displacement is given by 5Q{r) = SS{r)/TQ] we switch to Q simply because heat capacity traditionally is 
defined via this quantity. 

The isothermal bulk modulus (inverse isothermal compressibility) is defined by 



If di is the derivative with respect to the i'th spatial coordinate Xi [where {x,y,z) = {xi,X2, x^)] and Ui is the i'th 
component of the infinitesimal displacement vector u, the strain tensor e = is defined Q by 



diUj 



(20) 



The relative volume change is given by the trace of the strain tensor: 



SV 

_=tr(|) = V-u. (21) 

Denoting the stress tensor by a; = cr^ , the shear modulus G is defined [sl as follows (note that the isothermal and 
adiabatic shear moduli are always identical): 



(Ty - -tr(^)(5y = 2G (^cy - -tr(|)<5y j . (22) 

If there are both infinitesimal displacements and infinitesimal spatial temperature variations, which is generalized 
into the so-called Duhamel-Neumann relation [3, 0] that is the following constitutive relation linking mechanical and 
thermodynamic properties: 



<7y = KTtr{g)S,j + 2G i^e,j - -tr(|)% j - iSySTS,, . (23) 

For relaxing systems Kt and G are complex and frequency dependent, and for periodically varying boundary condi- 
tions both strain and stress tensors are generally complex. By the correspondence principle the Duhamel-Neumann 
relation applies also for the frequency-dependent case. 

The isothermal longitudinal modulus Mt is defined Q by 



Mt = Kt + -G. (24) 

Similarly, the adiabatic longitudinal modulus Ms is defined by Ms = Ks + (4/3)G where Ks = 1/ks is the adiabatic 
bulk modulus. As shown in detail in Ref. (23j (but implicit [36| already in Ref. Q), a useful quantity for the 
description of one-dimensional thermoviscoelasticity is the "longitudinal specific heat" defined by 



(25) 
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An important result of the present paper is that the longitudinal specific heat also plays a significant role for the 
spherically symmetric case. Again, c; is frequency dependent for relaxing systems. A useful identity psj follows from 
Eqs. (A.l) and (A.7): 

Ms Ks + lG 

From this one finds that c/ may be interpreted as a generally complex convex combination of Cp and cy '■ 

3Kt 4G 

This shows that the longitudinal heat capacity is effectively in between the isobaric and isochoric heat capacities. In 
analogy to the standard abbreviation 

7^^, (28) 
cv 

we define 7;, a quantity that is also generally complex and frequency dependent, as follows (23j : 

7/ = — • (29) 
Cv 

If the heat-current density is denoted by j, the heat conductivity A is defined via Fourier's law, 

j = -AV(5T . (30) 

The heat conductivity is generally assumed to be frequency independent, an assumption that was recently confirmed 
[33, [11]. The below solution, however, applies also if A were to depend on frequency. 

III. THE COUPLED EQUATIONS FOR THE TEMPERATURE AND DISPLACEMENT FIELDS 

A. The frequency-independent case 

As mentioned, entropy conservation always applies to first order, a fact that in Ref. Q is referred to as "the equation 
of continuity for heat." In terms of the infinitesimal displacement field u(r, t) and the infinitesimal temperature field 
(5T(r, i), for a volume element dV at average mass density po, the basic thermoviscoelastic equations of motion 
reflecting entropy conservation and Newton's second law, podVcPui/dt"^ — dVJ2j ^jcrij [2|, i3i, 1^, i23|] are (utihzing Eq. 

my- 

Po^ = MrViV • u) - GV X (V X u) - /SyV^T (31) 

and 

cv^+To(3v§-^iy-n) = XVHT. (32) 

We shall only be concerned with the inertia-free limit, i.e., the limit where the sample is much smaller than the 
wavelength of sound waves at the frequencies of interest. In this limit the acceleration term is negligible, thus 

MtV(V • u) - GV X (V X u) - Pv^ST = (33) 
cy^+ro/3y^(V.u) = XVHT. (34) 
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Before proceeding we note that the identity Pv = ctp Kt [Eq. (A5)] implies that if there is no thermal expansion 
upon heating, then /3y = 0. In this case the two equations decouple and reduce to the ordinary elastic equation of 
motion in the inertia-free limit and the heat-conduction equation, respectively. Thus the coupling between mechanics 
and thermodynamics arises only when the thermal expansion coefhcient is nonzero, as is indeed intuitively obvious. 
The pressure variation dp is defined by (assuming here and henceforth that the average stress tensor is zero) 

Sp^ -ltr{g). (35) 

This definition applies also for nonhydrostatic conditions, i.e., when a_ is not proportional to the unit tensor. When 
the specific heat or the expansion coefficient is termed isobaric, it refers to a situation where the trace of the stress 
tensor is constant, not necessarily its individual diagonal components. The relative volume change of Eq. (|2ip is 
related to changes in pressure, Sp, and temperature, 6T, by the following equation (see, e.g., Ref. [23| ) 

V ■u = --^Sp + ap6T . (36) 
Kt 

Using Eqs. (|36p and (Al) one finds that under isobaric conditions To/SyV • u = ToPyctpST — {cp — cy) ST. Thus Eq. 
(p4)) becomes the standard heat-diffusion equation 

^^DpV'ST, (37) 

where the heat-diffusion constant Dp involves the isobaric specific heat. Dp ~ X/cp. In general, isobaric conditions 
do not apply and the full coupled system Eqs. ([55)1 and must be solved. 

If V X u = 0, Eq. p3|) simplifies considerably. This applies for the spherically symmetric case, u — u{r)r, where 
Eq. ([55]) reduces to 

V(MtV • u - PvST) = . (38) 

This is immediately integrated to 

V .u^ ^6T + ai. (39) 

Here ai is an integration constant that is a function of time determined by the boundary conditions. Substituting 
Eq. ^ into Eq. ^ yields 

c,^+ToPy^^XV^6T. (40) 

B. The case of spherical symmetry and periodically varying fields 

We now assume that all fields depend only on radius r and that their time dependence is harmonic, i.e., oc exp(st) 
[s ~ ±iuj]. Henceforth u and ST refer to the complex frequency-dependent amplitudes of the infinitesimal radial 
displacement and temperature fields, respectively. All constitutive quantities are generally complex functions of the 
Laplace frequency s. 

With these assumptions, if differentiation with respect to r is denoted by a prime, using the identity V-u = r^'^{r'^u)' 
Eqs. ([391) and (gOl) become 

r'^{r^u)' ^^ST + a^ (41) 



and 
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cisST + TaPvsai = \r-'^{r'^5T')' . (42) 

Throughout the paper it is important to remember that not only the field amplitudes u and ST are generally complex 
functions of s, but so are all constitutive properties like (3vi Mt, etc. For simplicity of notation, however, the frequency 
dependence will usually not be explicitly indicated. 

We need two further equations for two auxiliary fields. One is the "radial pressure," 

5pr = —(Trr , (43) 

i.e., the normal force per unit area in the negative radial direction exerted on a spherical surface of radius r from 
the outside. The radial pressure is generally different from the pressure 5p of Eq. (|35p : the relation between the two 
follows from the definition of the shear modulus Eq. ((22|l : 

Sp^ ^Sp-2G (^err - iv ■ . (44) 

If G = - in particular for any liquid at zero frequency - radial pressure equals pressure. Since e^r = u' where e^r 
is the rr 'th component of the strain tensor in spherical coordinates Eq. (|44p and the Duhamel-Neumann relation 
Eq. (Ell) imply 

Sp^ ^-(Kt + -G)u' - (Kt - -G)-u + KraJT . (45) 
3 3 r 

The other auxiliary field to be introduced is the time-integrated heat-current density 5q. Invoking Fourier's law of 
heat conduction Eq. (PP]) we get 

6q ^-j = --6T' . (46) 
s s 

C. Scaling to dimensionless variables 

The solution of the thermoviscoelastic equations in spherical symmetry is extremely involved when expressed in 
terms of the original variables. We solve the equations below by proceeding in analogy to the solution of the one- 
dimensional problem detailed in Ref. [23^ where considerable simplification was obtained by scaling with suitable 
variables. One is free to scale by a characteristic length, time, temperature, and mass, or other combinations of these 
four fundamental dimensions. We choose pressure instead of mass as the fourth fundamental dimension. There is no 
mathematical requirement that scaling variables must be real; in the frequency domain they may well be complex 
and frequency dependent. 

The scaling units are chosen as follows. The unit of time is the inverse complex Laplace frequency, s^^. The unit 
of length is the (complex and frequency dependent) heat-diffusion length 

Id = \/^. (47) 
V QS 

It is occasionally convenient to use the (complex and frequency dependent) wave number k = 



The unit of temperature is the average temperature To, and the unit of pressure is the (complex and frequency 
dependent) isothermal bulk modulus K^- With these units we define the following dimensionless variables 
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i 


= st 


(49) 


f 


= t/Id 


(50) 


u 


EE u/Id 


(51) 


sf 


= ST/To 


(52) 


Spr 


= Spr/Kr 


(53) 


Sq 


= Sq/iKrlD) 


(54) 


c 


= TqCi/Kt 


(55) 


9 


= AG/{3Kt) 


(56) 


a 


= Toap . 


(57) 



It should once again be emphasized that all thermodynamic and mechanical constitutive quantities are generally 
complex and frequency dependent, although for simplicity of notation this fact is rarely explicitly indicated below. 

The scaled problem involves only the three independent frequency-dependent response functions: c, g, and a. 
Nevertheless, it is convenient to introduce two further dimensionless frequency-dependent parameters, 7 and 7; of 
Eqs. (p8|) and (|29p . which are functions of c, g, and a. Before proceeding we note that via the above definitions, Eqs. 
(A.l), (A. 4), and (A. 5), the parameters 7 and 7; obey 



I 52 

71 c(l -I- g) 

^ - 1 + -7?X,. (59) 

II 0(1+ g) 



These equations imply the following identity that turns out to be useful. 



7 + g = 7/ (1 + .9) ■ (60) 

When rewritten in terms of the above dimensionless variables, Eqs. (|42)l . ((4T|) . (|45| . and ((46)) become (where prime 
now implies differentiation with respect to f) 



r^'ir^STj = 5T+-ai, (61) 



c 



a 



f-2(f2u)' = -^^ST + ai, (62) 

6pr = -{l + g)u' -{2-g)f-^u + a5f, (63) 

5q = -dSf'. (64) 

These are the fundamental dimensionless thermoviscoelastic equations of spherical symmetry to be solved. 



IV. SOLUTION IN TERMS OF THE TRANSFER MATRIX 



Equation ([62]) leads to 



-^6f + ai^2f-^u, (65) 



by which Eq. (|63|) is simplified as follows 



(1 + g)ai + 3gf 



(66) 
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Equation implies (2/r)Sf' + Sf" = Sf + {a/c)ai, which is readily solved by substituting Sf = f{f)/r. The 
solution is 



ST— ai + f ^e^a^+f a^. 

c 



(67) 



where 03 and 04 are integration constants. Substituting this into Eq. (|62p and solving for u gives (where 02 is a 
further integration constant) 



1 / a 
u — - \ 1 — rr-r I r 



3V (1+5)5 



ai + f a2 + f (f ^-r 2)e^a3 - (f ^ + f ' ai] 
1 + a V / 



Now Eq. (|66p becomes 



Spr 



~ 2 

1 + ^ 1 ai + 35f-'^ a2 - ^-^ar'^ ((1 - f )e^ 03 + (1 + r) 



1+5 c 



1 + 5 



e 04 



Finally we have 



(68) 



(69) 



Sq = cf"2[(l - f)e'' 03 + (1 + f)e~'' 04] 
The solution is summarized in the form 



(70) 



Sf{r) 

svlf) 



as 
04^ 



(71) 



Here, if SQ{r) is the time-integrated heat current through the spherical surface of radius r, we have introduced 



SV 
SQ 



r^Sq 



SV 
4^ 



(72) 
(73) 



The matrix M(f) is given by 



/-(l + lfsf) ^~9f~' -3^af-3(l-f)e^ _3_1^5f-3(i + 



M(f) 



(1 - to)^" 



V 



'Tfl(l--V 
c(l — f )e'' 



Defining the transfer matrix by 



(74) 



T(f2,ri) = M(f2)M-i(fi), 
allows one to express the fields at in terms of those at f\ as follows 



(75) 



(Spr{f2)\ 

Sfih) 

sv{h) 



T(f2,fi) 



(Spr{fi)\ 

Sf{h) 



(76) 
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In this way the inner and outer boundary conditions of Fig. 1(b) are hnked by a matrix containing all relevant 
information about the physics of the system. The transfer matrix is 4 x 4, reflecting the fact that out of the eight 
variables of the four energy bonds of Fig. 1(b), four may be externally controlled. Via the transfer matrix the remaining 
four are determined by the four basic equations (|33p and (|34p . Once the transfer matrix has been calculated, various 
experimentally relevant special cases may be worked with modest efforts (Sec. V). 

By explicit calculation, e.g., via a computer program, one finds that T(f2,fi) has the form 



T(f2,fl) 



+ f-(f2,ri) 



-{f2-fl) 



+ T+(f2,ri)e^ 



(77) 



where 



I 7+gr 



1+a 

a 

c{i+g) 

3(7+9) 





7+9_ 

£(i+s) 

7+9 





£(1+9) 



a. 

£(1+9) 





(78) 



T-(f2,ri) 



/ 3a^g(l+r2)(fi-l) 
2c(l+g)2f^ 

a(l-fi) 
2c(l+g)f2 

2£(l+9)^ 
a(l+f2)(l-ri) 
2(1+9) 



3ag(l+r2)(l-fi) 
fl-1 

2f2 

a(l+f2)(l-ri) 

W+W) 
£(l+f2)(^i-l) 



9a^g''(l+f2)(l--ri) 
2c(l+g)2f|fa 
3ag(fi-l) 
2c{l+g)f2f'i 
3a''g(l+f2)(l-fi) 

3ag(l+f2)(fi-l) 
2(l+g)f? 



3Sg(l+f2) \ 
1 

2cf2ri 
a(l+f2) 
2c(l+g)fi 

l+f2 

2fi 



(79) 



and 



T+(f2,ri) 



' 3a^g(l-f2)(l+fi) 
2£(l+g)2fi 
a(l+fi) 
2c(l+g)f2 
a^(l-f2)(l+fl) 

2£(l+g)^ 
a(f2-l)(l+fl) 



V 2(l+g) 



3ag(r2 — l)(l+ri) 

2(i+g)fi 

1+ri 

2f2 

a(f2-i)(i+fi) 

2(i+g) 
£(l-f2)(l+ri) 
2 



9a 



'^g''(^2-i)(i+ri) 

2£(T+|p?pf 
3ag(l+fi) 
2c(l+g)f2?^ 
3a^g(f2-l)(l+fi) 

2£(T+Ip?f 
3ag(l-f2)(l+^l) 
2(l+g)r^ 



3ag(l-f2) \ 
2£(T+IFpr 
1 

2cf2ri 

5(1 — ^2) 

2c(l+g)fi 
r2 — 1 
2fi 



(80) 



Although the expression for the transfer matrix is cumbersome, T has a number of simple mathematical features. 
First, explicit calculation shows that det(M(f)) = 2c(l + g) independent of radius. This implies that 



det(T) = 1 



(81) 



Moreover, the definition of T immediately implies that T(f, f ) = E where E is the unit matrix, and that the inverse 
is given by 



T-l(f2,fi) =f (fi,f2). 

More generally, by its definition the transfer matrix clearly obeys the functional equation 



(82) 



T(f3,ri) =T(f3,?2)T(f2,ri). 
This gives rise to a number of relations between T", T~ and T+: 



(83) 



f°(r~3,fi) - f"(f3,f2)f"(f2,ri), 

f-(f3,ri) = f-(f3,?2)T-(f2,ri), (84) 

f+(f3,ri) = f+(f3,r2)T+(f2,fi), 



and 



f°(f3,r2)f+(f2,ri) = , T0(f3,?2)T-(f2,fi) = 0, 
f+(f3,?2)T°(f2,ri) -0 , T-(f3,f2)T0(f2,ri) = 0, 
f+(f3,?2)t-(f2,ri) -0 , T-(f3,f2)T+(r2,ri) =0. 

The dimcnsionless transfer matrix T(f2,ri) given by Eq. (|76|) is related to the dimensionless 
r(f2, fi) defined by the relation 



r(r2,fi) 

If T and F are split into four blocks of (2 x 2)-niatrices as follows 



/6pr{'fi)\ 
Sfih) 



( 5y{ri) \ 

SQih) 



and 



finds 



Calculating f from T gives 



Ti T2 
T3 T4 



ri T2 

f 3 f4 



-Tif7lf4 + T2 



r2 - 

f4 = 



-f 1T3-I 



rii(f2, fi) 

fl2(»'2,^l) 
fl3(»'2,^l) 
f22(?'2,^l) 



' 2 " ' 1 
Sail 

3(7 + g) 
^3 _ ^3 

3(7i - 1) 



ri 



f2 cosh(f2 — fi) — sinh(f2 — ^1) 



c(f| - fl) cfi[{f2 — fi) cosh(f2 - fi) + (fif2 — 1) sinh(f2 - ?i) 
3(7i-l) , 1 



c(f| - ff) c[(f2 - fi) cosh(f2 - fi) + (fif2 - 1) sinh(r2 - ^1) 
All off-diagonal elements of the four block matrices become identical, i.e., 



El 4 — Tti — r23 — — r34 — r4i — r43 — r 



Furthermore, 
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f3l(^2,?'l) = -f 13(?'l,?2) ^fl3{f2,ri), 

f 33(^2, ^^l) = -fll(7'l,?2) , 

f42(^2,f'l) ^ ~i'24{ri,f2) ^T24:{r2,ri) , 

f44(^2,f^l) = -f22(7^1,?'2) ■ 
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The response matrix F is symmetric, i.e., f i = T[, f 2 = fg and f4 = These Onsager relations follow also 
from general arguments (the fluctuation-dissipation theorem). Interestingly, T has an even higher symmetry: 

fl(f2,fl) = -f4(ri,f2) , f2(f2,ri) = -f3(fl,f2). (102) 

These relations follow from Eq. ([5^ and the connection between T and T. On the other hand, the symmetry relation 
r = r' implies 

T4T^ = T3T; , f^fi = T;T3 , E = f ;f 4 - , (103) 

and that 



TiT^ = T2T'i, T^T2=T'2T4, E = TiT^ - T2T;5 . (104) 

V. SOME CASES OF EXPERIMENTAL RELEVANCE 
A. A massive sphere 

In this first application (Fig. 2) of the formalism we inquire into how a solid sphere responds to a compression 
—SV and a heat supply —6Q [33| applied at radius ^2 (the transferred heat is positive when SQ < because 6Q refers 
to the heat flow in the positive radial direction). This is calculated from the transfer matrix by putting ri = 0. If 
*C \Id\, we expect the response matrix Go to be given by the constitutive equations (5), (7), and (2) generalized 
to frequency-dependent coefficients (where all variables refer to the outer radius r2): 



In general, the relation is 



Spr 
5f 



= G 



-SV 
-SQ 



(106) 



with 



G = -Ti(f2,0)T3-^(f2,0) = r4(f2,0) . (107) 

This follows from Eqs. ((76|) and (|87p via the boundary conditions SV{fi = 0) = and SQ{fi = 0) = 0. Returning to 
dimensional variables, Eqs. (|106p and p07p yields 



SpA _l(Ks ^\ f^SV 



(108) 



where Jd is the function of frequency defined by 



1/1 sinh(a;) 1 ,, , s. 



7( \3 a::cosh(a;) — sinh(a;) 

The frequency dependence of fu derives primarily from that oi Id- Note that fo^^ for w — > and fo ^ 00 for 
w — > 00. Asymptotically, one has fjy = r2/(S"/ilu) for a; — > 00. The function fjj describes how heat diffusion affects 
the measurement of the thermal and mechanical properties of a massive sphere that can only, of course, be accessed 
at the surface. 
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5Pr 



FIG. 2: Compression —SV and heat input ~SQ at the surface of a sphere give rise to changes Spr and ST in radial pressure and 
temperature, respectively, at the surface. For a finite sphere the response (Eq. (|108p ) is given by the constitutive properties 
Ks, as, and cy, as well as by heat diffusion. Eight kinds of input/output relations are discussed in the text. 

Below, several different thermoviscoelastic experiments on a massive sphere are considered. In principle there are 24 
such experiments [1^1, corresponding to the 24 coefficients discussed at the beginning of Sec. II: One may choose any 
of the four variables Spr, ST, ~dV, or —SQ at the outer radius r2 as the controlled input, any of the three remaining 
variables as the measured output, and any of the last two to be fixed. Since there are only three independent functions 
of frequency in the matrix of Eq. (|108[) . it is superfluous to discuss all these 24 experiments. It does make sense, 
however, to discuss more than just three cases, because the experimental challenges may vary from case to case. 

We discuss eight cases, corresponding in the low-frequency limit to the eight frequency-dependent thermodynamic 
response functions kt, us, cy, Cp,ap, as, Pv, and jSs- These cases are detailed below where, occasionally, relations 
from the Appendix are utilized. 

a. Compression with isothermal boundary: 



At low frequencies this converges to Kt/V. This result shows how one is limited upwards in frequency when attempting 
to do isothermal bulk modulus measurements. Because //) — > oo at large frequencies, the equation also describes the 
transition to the adiabatic bulk modulus above the characteristic heat-diffusion frequency lod defined by ujd = -D/r|, 
where D is the heat diffusion constant and ^2 the sample size. 

The bulk modulus can be measured in the frequency range, 1 Hz — 50 kHz by the so-called piezoelectric bulk modulus 
gauge (PBG) [lB|. This is a piezoelectric ceramic hollow sphere that may be filled with liquid. The ceramic shell 
has electrodes on the inside and outside and thus constitutes an electrical capacitor. Due to the piezoelectric effect 
the frequency-dependent capacitance - that can readily be measured - depends on the bulk modulus of the liquid 
which can be found after a calibration of the PBG. This device has a radius of 10 mm and since typical liquid heat 
diffusivities are of order 0.1 mm^/s, the characteristic heat diffusion frequency becomes 10~^s~^. Thus experiments 
performed with the PBG above 0.1 Hz can safely be said to be adiabatic, despite the fact that no special measures 
are taken to make the boundary conditions adiabatic. 

b. Adiabatic compression: 



Ironically, this boundary condition is difficult to achieve experimentally, whereas the isothermal experiment gives Ks 
at most frequencies. Thus Ks is easier to measure than Kt [l6j . 
c. Isochoric heating: 




(110) 




(111) 




(112) 



The low-frequency limit is l/{Vcv) giving the isochoric frequency-dependent specific heat. The high-frequency limit, 
however. 
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involves the longitudinal specific heat. Note that in this limit - even though the overall volume is constant - it is ci 
that appears, not cy. 

d. Heating into a free mechanical surface: 

=G..-%^ = -^fl + 7(/z.~l)). (114) 



-SQjp^ Gil Vcp 

The low-frequency limit is l/(Vcp) giving the isobaric frequency-dependent specific heat. The high-frequency limit, 

ST \ .,1 

for u! ^ oo , (115) 



-^Q J 47rr2YQA 

is identical to the isochoric high-frequency limit. Note that it is the longitudinal specific heat that enters into Eq. 
piSp . not Cp. This is similar to the fact that the frequency-dependent specific heat obtained from plane- wave effusivity 
measurements is not the isobaric specific heat, but the longitudinal '23l|. In that case nonisotropic stresses could be 
conceived as arising from the special kind of mechanical boundary conditions needed in order to keep the model of 
the plane-plate setup one-dimensional. Here we see, however, that nonisotropic stresses may arise in the liquid itself, 
not necessarily coming from clamping boundaries. This substantiates a conclusion of Ref. [2^, namely that it is 
not possible to probe the isobaric specific heat directly by effusivity measurements. - Note also that the thermal 
admittance per unit area is the same as for the planar geometry [23l | : 

Recently, radial heat effusion from the surface of a spherical cavity inside an infinite medium was shown also to involve 
the longitudinal specific heat f40]. 

e. Expansion by a controlled temperature oscillation at a free surface: 

i^^n .1 ^ = ^ (117) 



V G12G21 - G11G22 l + 7(,fo-l)' 
At low frequencies this approaches ap, whereas its high-frequency asymptotic form is given by 



^ '^^^ - ^ Y for Lu^ 00. (118) 



From Eq. p08)l it follows that 



V\STJ^^ (3ssr2 



V [STJ^^ " VTo \6pr 



(119) 



This is a radial version of the Onsager relation that corresponds to the Maxwell relation Eq. p^ . If an experiment is 
conceived where one measures the heat flux —6Q needed to keep temperature constant at the surface while applying 
a periodically varying radial pressure, one would find the same response function Eq. (jll7p . including the diffusion 
dependence. Similar radial versions of Onsager relations corresponding to the Maxwell relations Eqs. P7)) . and 

(fT7)) . hold for the last three response functions: 
/. Radial pressure in response to heating at constant volume: 



Spr \ 1 



SQjy VTaas 



(120) 



Thus as can be measured without interference from heat diffusion. This is not trivial, since the penetration depth 
\ljj\ of the temperature field into the sphere is frequency dependent, and it is the temperature field that creates the 
pressure variation. 
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g. Radial pressure in response to a controlled temperature oscillation at constant volume: This case leads to 



In contrast to case (f), this response function is diffusion influenced. It approaches /3y at low frequencies, whereas 



1 / 6pr 



- Y for Lo-.^. (122) 

V \ST J y assr2 



h. Volume expansion in response to heating for a free surface: 

SV \ 1 



SQJ„^ ToPs 



(123) 



As for case (f) we get a simple result that is independent of heat diffusion. 

In summary, the important role played by the longitudinal frequency-dependent specific heat is evident. Moreover, 
there is now an exact description of the transition from the adiabatic to the isothermal regimes of bulk modulus 
measurements utilizing the PBG [lB| ■ 



B. The "thermally thick limit" \Id \ <S r2 when ri <C r2 



We now proceed to discuss the case where both the inner radius is small and the sample is much larger than Ijj . This 
section prepares the theoretical basis of ongoing experiments where a small spherical thermistor placed in the center 
of the PBG makes it possible to simultaneously measure the frequency dependences of as , Ks and c/ on the same 
sample. Supplemented by shear modulus measurements [24] this provides a complete set of thermoviscoelastic response 
functions of a liquid. Below we determine the reduced transfer matrix X for a situation where a mechanical boundary 
condition at ri and a thermal boundary condition at r2 are stipulated. That is, X gives the linear relationship 



SQih)) -"^[sVir,)) ■ (124) 
There are four possibilities for X, denoted below by A, B, C, and D, depending on the different boundary conditions: 



X = A for SV{fi) = 0, SQih) = 0, 

X = B for 6V{h) = 0, Sf{f2) = , 

X = C for Sprih) = 0, <5Q(f2) = 0, 

X = D for 6prih) = 0, Sfih) = 0. 

If P is the inverse of T(r2, fi), i.e., P = T(fi, f2), one finds 



1 / ^21^32 — P22P31 P23P32 — P22P33 



^ I ^ zi^ ii i 2:2.^ ji ^ 2.3^ 32 22^ 33 I ('125) 

P32 \P41P32 - P42P31 P43P32 - P42P33J ' ^ ' 

1 / ^21^34 — P24P31 P23P34 — P24P33 



P34 \P4iP34 - P44P3I P43P34- P44P33J ' ^ ' 

C = — f^2lA2 - P22P1I P23P12 - P22Pl3\ (^27) 
P12 \P4iPi2 — P42P1I P43P12 — P42P13J 



£j _ 1 / ^21^14 — P24P11 P23P14 ~ P24P13 \ (128) 

Pl4 V^41 A4 - P44P1I P43P14 - P44P13 ' ' ^ ' 
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These expressions imply [4l| that 

det(A) = det(B) dct(C) = det(D) = 1 . (129) 

The exphcit expressions for A, B, C, and D are rather involved. We give only the components of the simplest one, 
A: 

^ ^ (r?+3(7i-l)(r;-fi)-r^)(r2 + l)exp(-(f2-ri)) + (r;-3(7i-l)(r;+^^ QgQ<> 

11 3afi7i((f 1-1) (7=2 + 1) exp(-(j=2-ri)) + (l-f2)(f 1 + 1) cxp(j=2-f i)) ' V ) 



A12 

Ml 
A22 



(g(f? + (7i-l)(3f?-3fi+r|))+7,f|)(f2 + l)cxp(-(f2-fi)) + (§(f;-(7,-l)(3f;+3fi-fi))+7,fi 
r^ari7i((fi-l)(f2+l) oxp(-(f2-f i)) + (l-r2)(ri + l) oxp(f2-ri)) 

3a7! ' 
c(g''i+7f^) 



The transfer matrix A becomes much simpler when ri <^ r2 and \Id\ <C r'2. In terms of the scaled variables we seek 
the limits ri/f2 ^0 and \ f2\ — > 00. In these limits one finds from Eq. (|130p that A ^ © where 
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_ I 37iQfi(l+fi) 7!afi(l+ri) j flSll 

\ 3Q7i Q7; / 

By explicit calculation one finds that B, C, as well as D, all converge to in the same limits. That is, in these limits 
the linear connection between the thermal response at the inner radius and the mechanical stimulus at the outer 
radius is independent of the mechanical boundary condition at the inner radius or the thermal boundary condition at 
the outer radius. If one inverts A, B, C, and D and go to the same limits, the inverse matrices all converge to the 
matrix ^ given by 



* = [^h ■ (132) 

\ 3Q7i 3a7ifi(l+fi) / 

In order to make the above results more transparent we return to dimensional variables. Thus introducing the 
sphere volume 

Air 

V2 = -jvl, (133) 
and the quantity Zth to be identified below with a thermal impedance, 

Ztx. = 7 ^ , (134) 



47rAri (l + ^srici/X^ ' 



Eq. (fmj) becomes 



and Eq. ()132|) becomes 



(5T{ri)\ _ fsZthV2Toas sZthToPs\ fSpr{r2)\ /-,oc\ 
[6Q{n)) - [ V2Toas To/3s ) [sv{r2)j ' ^'"^^^ 

f5pr{r2)\ _ ( TqPs -sZthTo(3s\ fST{ri)\ , . 

\SV{r2)J \-~V2T0as sZthV2ToasJ \SQ{n)J ■ ^^''''^ 

Three results may be inferred: 1) The adiabatic frequency-dependent bulk modulus is given by the normal-stress 
response to a compression at the outer surface independent of the other three boundary conditions: 
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2) The temperature response on the inner surface to a thermal current s5Q is the thermal impedance involving the 
longitudinal specific heat. This again holds independent of the other three boundary conditions: 

^)(.r) = Z,,. (138) 

3) The relation between the heat displacement at radius ri and the volume and negative normal-stress variations at 
radius r2 does not involve any "delay" caused by thermal diffusion: 

5Q{ri) = ToV2as Spr{r2) + ToPs SV{r2) . (139) 

C. Mechanical boundary conditions 

If we control the boundary conditions solely via isobaric or isochoric constraints, there is a thermal transfer matrix 
Y connecting heat and temperature variations at the inner and outer radia: 

^^(-)^-Y(...ors^)V (140) 



(141) 
(142) 
(143) 
(144) 



Depending on boundary condition we define [4 







H(f2,ri) 


for 


SVih) 


= and SV{f2) = 


0, 


Yih: 


^h) = 




for 


SV{h) 


= and 5pr(f2) — 


0, 


Y(f2, 


^h) - 


L(r2,fi) 


for 


Spr{fi) 


= and 5V{f2) = 


0, 






N(f2,ri) 


for 


Spr{fi) 


= and 5pr{f2) = 


: 0. 


all four cases 


one has 











det(Y) = 1 . (145) 



The following identities follow trivially 



H-i(f2,ri) = H(ri,f2) (146) 
K-i(f2,ri) = L(fi,f2) (147) 
N-i(f2,ri) - N(ri,f2). (148) 

Motivated by an experimental setup studying a liquid drop on a thermistor bead [20j we calculate K explicitly (see 
Figs. Oandd]). In analogy to the calculation leading to Eq. (|125p we find in terms of the matrix elements of T (where 
all T matrix elements are evaluated at (f2,fi)) 



K(f2 ri) — ^ f T22T11 — T21T12 T24T11 — T21T14 
Til V^42T'ii — TiiTi2 Ti4:Tii — TiiTi4 



(149) 



Switching to real units, K is defined via 



mr2)\ _ ^ fSTin 
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FIG. 3: AC-calorimetry in spherical geometry [23|. A liquid droplet of of diameter 0.7 mm is placed onto a strongly temperature- 
dependent resistor of radius 0.3mm. By the so-called 3a; technique the apparent heat capacity Capp = 5Q\/ST\ was found. At 
low frequencies in a narrow frequency range this gives to a good approximation the isobaric specific heat via Eq. (|157p . At 
higher frequencies heat diffusion makes the extraction of the frequency-dependent specific heat more involved. 



5Vi=0 5Pr,2 = 
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5Ti 5Qi 5T2 5Q2 



Thermistor Liquid Surroundings 

FIG. 4: Energy bond graph diagram 'ss'] of the physical interactions between the thermistor, liquid, and thermal leak to the 
cryostat, modeling the setup of Fig. 3. The thermistor acts as both heat source and thermometer. The liquid is described by 
the transfer matrix K for the case with mechanical clamping at the thermistor and a free outer surface. 



and related to K by (where k = l/^u) 



ifii(r2,ri) = kii{kr2,kri) , (151) 

^I2(r2,ri) = ki2{kr2Mi) , (152) 

K2i{r2,n) -^K2i{kr2,kn), (153) 

K22{r2,ri) = k22ikr2,kn). (154) 

We shall not exphcitly give the components for the general case, but limit ourselves to the thermally thin limit 
r2 ^ \Id\ (i.e., \kr2\ ^ 1) where the results simplify considerably. In this limit one finds 20] after Taylor expanding 
the matrix elements of T that 



where, if y = (47r/3)(r| — rf) is the volume, 



R 



1 



1 



1 



AttX \ri r2 



(156) 



and 
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ifT(r2/ri)3 + |G 



C = cvV , (157) 



play the role of thermal resistance and capacitance, respectively. These results imply that in the thermally thin limit 
Cp is measured if ^2 ^ ri, whereas ci is measured if r2 = ri. We finally note that the above formulation is easily 
incorporated into a model taking into account the thermal heat loss to the surroundings [20j . 



D. No thermomechanical coupling 

As mentioned in the beginning of Sec. Ill, if the isobaric thermal expansion coefficient ap is zero, one has j3v = 
and there is no thermomechanical coupling. In this case, the components T21, r23, T41, and T^^i vanish, implying that 
heat and temperature variations at f2 only depend on heat and temperature variations at fi. Also, T12, T'14, T32, and 
T34 vanish, implying that pressure and volume variations at f2 depend only on pressure and temperature variations 
at fi. Note that when there is no thermomechanical coupling, all specific heats are identical: 

Cp = cv = ci . (158) 

This is often a good approximation for solids, but rarely for liquids 

When there is no th 
defined [i^ as follows. 



When there is no thermomechanical coupling, heat diffusion is described by a 2 x 2 thermal transfer matrix T*^ 



The components of T**^ are found by substituting Op = into Eq. ([77)1 . The results arc as follows (where k = I/Id) 

Tl^ = — cosh (/c(r2-ri)) + -^ sinh (fc(r2-ri)), (160) 
r2 kr2 

Tl^2 = (161) 



^21 



47rA krir2 



(1 — k^rir2) sinh {k{r2 — rij) — k{r2 ~ r\) cosh (fc(r2 — ri)) , (162) 

T^^ = cosh (/c(r2 -ri)) - sinh (fc(r2 -r-i)) . (163) 
r\ kr\ 

(164) 

Interestingly, the same purely thermal 2x2 transfer matrix describes a low-viscosity liquid (g 0) even when ap ^ 0, 
if either the inner or the outer surface is free, i.e., if 5pr{ri) = or if 5pr{i'2) = 0. In these two cases, however, not all 
three specific heats are identical, only Cp = c; applies. This is consistent with the above remark regarding the validity 
of the standard heat diffusion equation ((37|) . 



VI. CONCLUDING REMARKS 



Thermoviscoelastic response functions are notoriously difficult to measure. This paper establishes the theoretical 
framework necessary for developing experimental methods that utilize spherical symmetry for measuring such response 
functions. From the complete solution of the problem in the form of the transfer matrix the equations describing any 
realistic experimental situation may be derived, as exemplified in the last section. 

The thermoviscoelastic response functions are important to determine for liquids approaching the glass transition 
(still in metastable equilibrium above the transition). For such ultraviscous liquids all thermodynamic coefficients 
become complex and frequency dependent for frequencies in the range of the inverse Maxwell relaxation time. To 
the best of our knowledge there are yet no reliable measurements of a complete set (i.e., three [11, 0, IH, [28j ) 
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of thermoviscoelastic response functions for any such liquid. The determination of such complete sets, from which 
all other thermoviscoelastic response functions are easily calculated, serves the obvious purpose of elucidating the 
macroscopic dynamics and thermodynamics of ultraviscous liquids. Very recent theoretical developments even further 
stress the importance of developing reliable methods for measuring thermoviscoelastic response functions. It now 
appears that the class of van der Waals liquids (possibly supplemented by some liquids forming bulk metallic glasses) 
have particularly simple properties: Liquids with non-directional chemical bonds exhibit strong correlations between 
equilibrium pressure and energy fluctuations (28l . |4^ . For such "strongly correlating viscous liquids" it has been 
shown that there is basically only one independent thermoviscoelastic response function 25] . It would be interesting 
to have this prediction subjected to experimental tests. Moreover, for strongly correlating viscous liquids there are 
strong indications from computer simulations that thermoviscoelastic measurements can determine the exponent of 
the so-called density scaling that collapses the relaxation time's pressure and temperature dependence onto a master 
curve [isj . 

As regards the above results, it is notable that the longitudinal specific heat c; plays a dominant role. It is perhaps 
not surprising that c; enters repeatedly into the equations describing the one-dimensional case [l^] - after all, this is 
what it was defined to do - but it is less obvious that c; also plays a dominant role for the case of spherical symmetry. 
Note, however, that this result is implicit already in Ref. [3|]. 
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APPENDIX 

If (^p = Cp/To and = cy /Tq, several standard thermodynamic relations between the eight linear-response functions 
CpT Cvt i^T, Ks, ctp, as, Pvy Ps are summarized below. Via the correspondence principle all relations apply also for the 
corresponding complex, frequency-dependent linear-response properties. Note that the inverse of the isothermal bulk 
modulus Kt is the compressibility, l/Kx = kt, and similarly for the adiabatic bulk modulus, 1/Ks = ks ■ 



ktPv = appv (A.l) 
' ' (A.2) 
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